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O ' Abstract 
• 1— I " 

Analytically and quantitatively we reveal that the GLE equation, based on a mem- 
ory function approach, in which memory functions and information measures of 
statistical memory play fundamental role in determining the thin details of the 
stochastic behavior of seismic systems, naturally conduces to a description of seis- 
J> ' mic phenomena in terms of strong and weak memory. Due to a discreteness of seismic 

■ signals we use a finite - discrete form of GLE. Here we studied some cases of seis- 
j- , mic activities of Earth ground motion in Turkey with consideration of complexity, 

CO \ nonergodicity and fractality of seismic signals. 
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1 Introduction 



Specific stochastic dynamics occur in a large variety of systems, such as su- 
percooled liquids, seismic systems, human brain, finance, meteorology and 
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granular matter. These systems are characterized by an extremely rapid in- 
crease or a slowdown of relaxation times and by a non-exponential decay of 
time-dependent correlation functions [ip]. 

The canonical theoretical framework for stochastic dynamics of complex sys- 
tems is the time-dependent generalized Langevin equation (GLE) [3llHl5]|6irflll4fl5l 
It successfully describes the phenomenon of statistical memory, whereby the 
relaxation time for order parameter fluctuations scales as a power of the corre- 
lation length. An obvious question to ask would be whether this framework can 
be adapted to describe seismic phenomena. Analytically and quantitatively we 
show that the GLE equation, based on a memory function approach, where 
the memory functions and information measures of statistical memory play 
fundamental role in determining the thin details of the stochastic behavior of 
seismic systems, naturally leads to a description of seismic phenomena in a 
terms of a strong and weak memory. Due the discreteness of a seismic signals 
we use a finite - discrete form of GLE. Here we study some cases of seismic 
activities of Earth ground motion in last years in Turkey with consideration 
of complexity, irregularity and metastability of seismic signals. 



2 Some extraction from the theory of discrete stochastic processes 

The GLE analytical model was originally proposed for displaying stochastic 
behavior of signals in complex systems [3llHl5]|6ll7] . one of which identifies mem- 
ory effects with diverse memory time scales in manyfold of signal's correlation 
so that the arbitrary seismic state is characterized by the set of memory time 
length scale in a system. 

Here we consider data of seismic signals recording as a time series ^: 

^ = {^0, ei, 6, • • • , ^N-i} = {^(0), e(r), e(2r), . . . , ^([iV - l]r)}. (1) 

Here r is a discretization time of seismic signals, is a total number of signals. 
A set of fluctuations 6C, is an initial dynamic variable Wq : 

1 Af-l 

Wo = {6^0,5^1, 5^2,..., 5^N-i}, 50 = 0-(0, (0 = T7E^r (2) 
The Gram-Schmidt orthogonalization procedure 

{Wn,W^) = 6n,m{\Wn\^) (3) 
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leads to the set of the following orthogonal dynamic variables: 













W2 




- AiW^o, 



_ Wn+l = CWn - KWn-1, n>l, 

where £ = (A — 1)/t is the Liouville's quasioperator and A„ is the relaxation 
parameter of the nth order (where A is the shift operator Axj — Xj^i and r 
is the discretization time). 

Within the framework of statistical theory and Zwanzig-Mori's theoretical- 
functional procedure of projection operators one can receive following recur- 
rent relation as a finite-difference kinetic equation: 



m—l 

AM^{t) = TA„+iM„(i) - T^A^+i Yl ^n+i{t - jr)M^{jT), 0,1,2,.... 

j=0 

(5) 

Here we introduce a Liouville's quasioperator eigenvalue A„+i, a relaxation 
parameter A„+i and a memory function M„(i) of the nth order respectively 



An 



Ar,. = 



(l^-n M (t) = (6) 



For analysis of relaxation time scales of underlying processes we use the 
frequency-dependent statistical non-Markovity parameter e„(a;): 



;^n-l(^) 



1/2 



(7) 



Here /ini^) is a frequency power spectrum for the memory function of the nth 
order: 

/x„(a;) = tJ2 Mnijr) cos{jTu) . (8) 

j=0 



Using of Eqns. (1) - (8) we can study all specific singularities of the statistical 
memory effects in an underlying system. Non-Markovity parameter and its 
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Fig. 1. An initial row data of seismic signals from 7 areas of seismic signals: grsn, 
kelt, mack, sgkt, uldt, seyt, gdz. The discretization time is r = 0.02 sec. 

statistical spectrum were introduced by Yulmetyev et al. in [8j. It is worthy 
of mentioning that non-Markovian character of seismic data was discussed by 
Varotsos et ah [H]. One of the first proofs of non-Markovity of empirical random 
processes was given in Refs. [TU]. Stochastic origins of the long-range correla- 
tions of ionic current fluctuations in membrane channels with non-Markovian 
behavior were studied in [TT] . 



3 An analysis of results 



Fig. [T] presents the initial time series of seismic signals for 7 seismic origins: 
grsn, kelt, mack, sgkt, uldt, seyt, gdz. Discretization time is r = 0.02 sec. 
From the Figures we can see that all time series have distinctive features. 

Fig.[2]demonstrates the frequency dependence of the first point of non-Markovity 
parameter for 7 seismic origins from Turkey £1(0;): grsn, kelt, mack, sgkt, uldt, 
seyt, gdz . Since the nature of each seismic source is unknown to us, it would 
be interesting to establish its character. It seems possible that the signals can 
be distributed to 3 groups: group A {kelt, gdz), group B {grsn, sgkt, uldt, 
seyt) and group C {mack). 
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Fig. 2. The frequency dependence of the first point of non-Markovity parameter 
for the each seismic origin: grsn, kelt, mack, sgkt, uldt, seyt, gdz. 

Signals of Group A are characterized by the more regular structure and 
smooth decay of the function ei{uj). 

There is a frequency dependence of non-Markovity parameter for seismic ori- 
gins: gdz kelt in Fig. [31 General power dependence Si^u) = {ijj/ujq)~°' is 
submitted by a continuous line with parameters ujq = 0.2, a = 1.05. 

Signals of Group B are characterized by the irregular frequency structure 
and by the frequency bursts on the distinct frequencies. Spectra have a noisy 
character. 

There is a frequency dependence of ei{u) for 4 seismic origins: grsn, sgkt, uldt, 
seyt in double log-log scale has been presented in Fig. 4. Discretization time 
is r = 0.02 sec. General power dependence ei{uj) = (cj/cuq)"" is submitted by 
a continuous line with parameters uq = 0.2, a = 0.6. 

Signals of group C can not be attributed to one of the above groups. The 
parameter ei{uj) fluctuate strongly between 1 and 10 in full frequency scale. 
That testifles existence of strong memory effects in the long-range time correla- 
tion. A possible origin of the similar signals is due to strong Earthquake. More 
careful and detailed analysis of the signal structure on the various time scales 
and relaxation levels (with the taking into account of the long-range correla- 
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Fig. 3. (Color online) The frequency dependence oiei{uj) in double log-log scale for 2 
seismic origins: gdz kelt. The time discretization is r = 0.02 sec. The general power 
dependence of = (lo/ljo)"" is submitted by a continuous line with parameters 

= 0.2, a = 1.05. 




Fig. 4. (Color online) The frequency dependence of £i{io) for 4 seismic origins: grsn, 
sgkt, uldt, seyt in double log-log scale. The time discretization is r = 0.02 sec. The 
general power dependence of £i{uj) = {u/loo)~°' is submitted by a continuous line 
with parameters ujq = 0.2, a = 0.6. 
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Fig. 5. (Color online) The frequency dependence of non-Markovity parameter £i{uj) 
for seismic origin mack in the double log-log scale. The time discretization is 
r = 0.02 sec. The general power dependence of eo{uj) = (w/wq)"" is submitted 
by a continuous line with parameters loq = 0.2, a = 0.4. 

tion, memory effects, nonergodicity and metastability of underlying system) 
is required. 

Fig. 5 displays the frequency dependence of ei{u) in double log- log scale for 
mack. Discretization time is r = 0.02 sec. General power dependence ei(uj) = 
[ui /ojq)~'^ is submitted by a continuous line with parameters ujq = 0.2, a = 0.4. 

The analysis of all spectra shows that all the signals can be classified into 
three different groups in the order of the breaking of fractal behavior of high 
frequency dependence of £1(0;). Signals for the group A can be characterized 
by stronger fractality with the exponents a = 1.05. A linear trend with the 
small fluctuation has been simultaneously observed in the spectrum. We can 
see range of the diversity 10 < Bi{uj) < 100. Signals for the group B are 
characterized by the breaking of fractality with the exponents a = 0.6 . A 
nonlinear oscillating trend with big fluctuation has been observed here. Spectra 
of signals for group C are characterized by a weak fractality with the exponent 
a = 0.4 and ei{uj) ~ 1. 



The auto-correlation functions (ACF) 
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for the signals of different groups have been presented in Fig. 6. 
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Fig. 6. (Color online) The auto-correlation function of the signals. 

The left panel includes ACF for the signals of group A and group C, {mack) 
and {kelt, gdz), while the right panel contains ACF for the signals of group B 
{grsn, sgkt, uldt, seyt). The brackets here note averaging in time iterations. It 
is seen from the Figures, that auto-correlation of the signals has a pronounced 
nonergodic character (undamped behavior of the time correlation function at 
time t — > oo): 



lim C{t) ^ 0. (10) 

t— >oo 

According to recent works on the ergodic hypothesis Eqns. (9), (10) would 
imply violation of ergodicity. Net results [12] suggest the breaking of ergod- 
icity for a class of generalized, Brownian motion, obeying a non-Markovian 
dynamics being driven by a generalized Langevin equation. This very feature 
originates from vanishing of the effective friction. Khinchins theorem of er- 
godicity is examined [T3] by means of linear response theory. The resulting 
ergodic condition shows that, contrary to the theorem, irreversibility is not a 
sufficient condition for ergodicity. 

Similar behavior of the time correlation functions is characteristic for the 
supercooled and glass states of condensed matter [14fl5]. A higher level of 
nonergodicity corresponds with the signals of group A and group C, whereas 
signals of group B are characterized by the minor nonergodicity. Signals from 



8 



Table 1 

The frequency relaxation parameters for seismic signals. 



Object 


Ai 


A2 


As 


Ai 


A2 


grsn 


-0.0306 


-0.2004 


-0.8945 


0.0497 


0.0396 


kelt 


-0.0667 


-1.1399 


-1.0485 


-0.0193 


0.2954 


mack 


-0.3938 


-0.5808 


-1.0747 


0.4374 


-0.0342 


sgkt 


-0.0306 


-0.7490 


-1.0849 


0.0156 


-0.2967 


uldt 


-0.0349 


-0.2605 


-0.7609 


0.0526 


0.1243 


seyt 


-0.0745 


-0.2434 


-0.8427 


0.1173 


0.0590 


gdz 


-0.0646 


-0.8836 


-0.9650 


0.0155 


0.2586 



the object uldt are rigorously ergodic. Therefore one can suppose that these 
signals cannot belong to Earthquake. 

Table 1 presents a set of relaxation parameters Ai,A2,A3,Ai and A2 for seis- 
mic signals from: grsn, kelt, mack, sgkt, uldt, seyt, gdz. The set of these 
parameters characterizes some peculiarities of relaxation processes on the low 
relaxation levels of seismic systems (1, 2 and 3). It has seen from the Table 
that these parameters don't have a clear distinction vs distinctions unlike those 
visible from the frequency dependence of ei{uj). Therefore thinner and more 
sensitive techniques are needed for studying of dynamic processes in examined 
signals. 

Figure 7 depicts the comparison of seismic data with results of computer 
simulation for the metallic glass AI<^qCu<^q. On top, auto-correlation functions 
of the seismic events {grsn, sgkt, uldt and seyt) are shown. At the bottom, 
auto-correlation functions for incoherent scattering of copper atoms are shown 
for comparison. The data are received with the help of statistical averaging on 
ensembles of statistical systems. Each curve was received by time averaging. 
Full sample consists of 45000 points. The single calculation was carried out 
for each separate time window of 500 points size with time step of 0.02 s. 
Further this window was displaced to one step to the right up to the end of 
time sampling. 

On comparison data for seismic phenomena with the results of computer sim- 
ulations for the glassy system it is visible, that in behaviour of correlation 
functions for EQ's nonergodic effects, characteristic for glass-hke behaviour of 
dense systems are distinctly observed. 

To illustrate the general picture of nonergodic singularities in chaotic seis- 
mic systems the parameter of nonergodicity / = lim c(t) for a set of seismic 

t— >oo 

events has been calculated. Parameter of nonergodicity for a set of seismic 
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Fig. 7. (Color online) On top, auto-correlation functions of the seismic events grsn, 
sgkt, uldt and seyt in time log-scale. At the bottom: temperature dependence of 
incoherent scattering function Fs{k,t) for the Cu - component in metallic glass sys- 
tem ^/soCuso-alloy for the wave vector k = 3.05A at the temperatures: T=2000 
K, 1000 K, 600 K, 500 K, 400 K, 200 K (bottom up). 



events appeared equal: 0.9995 (gdz), 0.9995 (kelt), 0.99885 (mack), 0.9310 
(grsn), 0.7241 (sgkt), 0.3965 (uldt), 0.0603 (seyt). The resulted data are ev- 
idence of strong singularity of seismic phenomena in 5 sources (gdz), (kelt), 
(mack), (grsn), (sgkt), they show moderate nonergodicity for a source (uldt) 
and weak singularity for a source (seyt) . All taken together received data well 
speak about wide variety of effects of nonergodicity in the seismic phenomena. 
Similar variety of nonergodicity effects can be very useful and extremely effec- 
tive for the classification of wide variety of seismic phenomena. Note, that the 
notions of fractality and non-Markovity have quite extensively been explored 
in the past, see, for example, [TB|17fl8|[TU] . 
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4 Summary 



In this work we have presented the results of statistical analysis of seismic 
signals in Turkey for 7 objects {grsn, kelt, mack, sgkt, uldt, seyt, gdz). Our 
study was made in the context of statistical theory on the discrete non-Markov 
processes, which is based on the Generalized Langevin Equation (GLE). It al- 
lows to take into account the effects of the statistical memory, metastability 
and space-time nonlocality. We have shown with the theory that all consid- 
ered signals can be divided into three groups in order of breaking of the frac- 
tal behavior in high frequency zone of spectrum of non-Markovity parameter. 
Signals from A group {kelt, gdz) can be characterized by the pronounced frac- 
tality, signals from the B group {grsn, sgkt, uldt, seyt) can be characterized 
by the moderate fractality and signals from the C group {mack) correspond 
to weak fractality and powerful non-Markov processes. From the analysis of 
the time correlation function we can confidently certify hypothesis Abe pQ 
of nonergodic "glass-like nature" of seismic signals for Earth activity. On the 
other hand aforementioned testifies an information concerning a wide variety 
of metastability in seismic phenomena. 
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